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We revise the symmetry analysis of a modified system of one-dimensional shallow-water equa- 
tions (MSWE) recently considered by Raja Sekhar and Sharma [Commun. Nonlinear Sci. Numer. 
Simulat. 20 (2012) 630-636]. Only a finite dimensional subalgebra of the maximal Lie invariance 
algebra of the MSWE, which in fact is infinite dimensional, was found in the aforementioned paper. 
The MSWE can be linearized using a hodograph transformation. An optimal list of inequivalent 
| one-dimensional subalgebras of the maximal Lie invariance algebra is constructed and used for Lie 

■ reductions. Non-Lie solutions arc found from solutions of the linearized MSWE. 

m 

■ 1 Introduction 

■ The computation of exact solutions of nonlinear systems of partial differential equations remains 
an important task, despite the increased interest in the numerical simulations of such differential 

,^Zh \ equations. Exact solutions are crucial, for instance, in the testing of numerical methods as well 

as in the analytical study of the associated partial differential equations. Finding exact solutions 
remains rather challenging, unless the equations of interest belong to a category for which well- 
known solving methods exist. However, in general, the most one can do is to construct particular 
solutions. Particular solutions, albeit often simple, can be of great help to gain understanding 
!> ■ and insight into the principal dynamics of such systems of equations. 

One way to obtain exact solutions of differential equations is through the study of their Lie 
symmetries. It is worth to know that there is also a number of non-Lie methods for finding 

■ exact solutions, see e.g. [10,15,24]. Lie symmetries and the related computational methods for 
. finding group- invariant solutions of systems of differential equations constitute a well investigated 

subject. The current paper concerns itself with a hydro dynamical problem. Relevant results on 
symmetries and exact solutions for hydrodynamical systems can be found e.g. in the textbooks [1, 
2,9,20,24,26,27] and in the papers [3,11,14,16,17,22,23,29-31]. 

The recent paper [34] intended to find Lie symmetries and exact solutions of the following 
system of modified one-dimensional shallow-water equations: 
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Ai := ut + uu x + g (l + ^ J h x = 0, A 2 := h t + uh x + hu x = 0, (1) 

where u is the fluid velocity in x-direction, h is the height of the water column, and g and H 
are constants related to the gravity acceleration and momentum transport, respectively. The 
constant g can be scaled to g = 1 by an equivalence transformation. In the case of H = 0, this 
system reduces to the usual form of the one-dimensional shallow-water equations. 

We found a number of incorrect and incomplete results in [34], which is why we reconsider 
the symmetry analysis of the system (1). One main concern is that the maximal Lie invariance 
algebra of (1) is supposed to be infinite dimensional. At the same time, in [34] it was found to 
be finite dimensional. This sits ill with the fact that system (1) can be linearized via a hodo- 
graph transformation and thus should display an infinite dimensional maximal Lie invariance 
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algebra. Another problem with the Lie reductions carried out in (1) is that no optimal list of 
inequivalent subalgebras was determined which should be the base for an efficient computation 
of group-invariant solutions. As a consequence, the exact solutions found in [34] are in an overly 
complicated form. Moreover, we explicitly show that one of the reductions in [34] is incorrect. 

The structure of this paper is as follows. In Section 2 we compute the Lie symmetries admit- 
ted by system (1). We find the infinite dimensional maximal Lie invariance algebra of system (1) 
and compare it to the algebra found in [34] which is only finite dimensional. We indicate how the 
modified shallow-water equations are linearized through a hodograph transformation. In Sec- 
tion 3 we classify inequivalent one-dimensional subalgebras of the infinite dimensional maximal 
Lie invariance algebra of the modified shallow-water equations. One-dimensional subalgebras 
of the symmetry algebra found in [34] are also classified. In Section 4 we construct the group- 
invariant solutions based on the optimal list found in Section 3. We also determine non-Lie 
solutions of the modified shallow-water equations through solutions of the linearized modified 
shallow-water equations. The results obtained are summed up in the final Section 5. 



2 Lie symmetries of the modified shallow-water equations 

In this section, we compute the maximal Lie invariance algebra q of the system (1). This is done 
by determining the coefficients of the infinitesimal generator 

Q = r(t, x, u, h)d t + x, u, h)d x + r](t, x, u, h)d u + <p(t, x, u, h)dh, 

of a general one-parameter Lie symmetry group. These coefficients are found using the infinites- 
imal invariance criterion, which in the present case reads 

Q«Ai = 0, Q«A 2 = 0, (2) 

where this equality has to hold on the manifold defined by Ai = and A 2 = 0. In system (2), 
denotes the first prolongation of the operator Q, which is of the form 

Q(i) = Q + rfd ut + rfd Ux + fid,* + <P x d hx . 

Here the coefficients are given by 

77* = D t n - u t D t T - u x D t £, if = D x n - u t B x T - u x D x ^, 
4> l = B t 4> - h t D t r - h x D t £, (f = B x <t> - h t B x r - h x B x t 

which follow from the standard prolongation formula, see e.g. [2,9,10,26,27]. The total derivative 
operators and D x arising in the expressions of the coefficients if , n x , (ft and <fi x of the prolonged 
operator are given by 

D t = d t + u t d u + h t d h + u tt d Ut + h tt d ht + u tx d Ux + h tx d hx H , 

B x = d x + u x d u + h x d h + u tx d ut + h tx d ht + u xx 0u x ~\~ h xx d/ lx -)-■■■. 

Explicitly, the infinitesimal invariance criterion (2) reads 

rf + urf + nu x + (l + j\ <jf - ^(f>h x = 0, 

4? + ucp x + nh x + hrf + <j)u x = 0. 

Plugging the coefficients (3) into the above system, substituting u t = —(uu x + (1 + H/h)h x ) 
and ht = —{uh x + hu x ) wherever they arise and splitting the resulting equations with respect 
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to powers of the derivatives of u and h, we derive the following system of determining equations 
for coefficients of the vector field Q: 




(4) 



4>t + u(j> x + hr] x = 0. 

This system has two inequivalent cases for solutions, depending on whether H ^ or H = 0. 

If H = 0, system (1) reduces to the usual one-dimensional shallow- water equations for which 
the general solution of the determining equations (4) is 



where c±, . . . , C4 are arbitrary reals constants and the functions / and g run through the set of 
solutions of the system 

9u ~ uf u + hfh = 0, g h - uf h + f u = 0. 

This is the linearized system of one-dimensional shallow-water equations, obtained from the 
shallow- water equations through the hodograph transformation in which f = t and g = x serve 
as the new unknown functions and u and h become the new independent variables. The maximal 
Lie invariance algebra of the one-dimensional shallow-water equations is thus spanned by the 
vector fields 

V 1 =td t + xd x , G = td x + d u , V 2 = 2hd h + ud u - td t , 
C = Ahud h + (4hg + u 2 )d u + (— 6ttt + 2x)d t + (6hgt - 3u 2 t)d x , 
£(f,g) = f(u,h)d t + g(u,h)d x , 

see also [7, 19]. 

We will only concentrate on the case H 7^ subsequently, which corresponds to the modified 
shallow-water equations. If H 7^ 0, then the system of determining equations (4) yields the 
general solution 

r = c\t + f(u,h), £ = c\x + c 2 t + g(u,h), rj = c 2 , <f> = 0, 

where c\ and C2 are arbitrary real constants and the functions / and g run through the set of 
solutions of the system 



r = c\t + Ci{2x — 6tu) + f(u, h), £ = c\x + c 2 t + C3X + C4t(6h — 3u 2 ) + g(u, h), 
rj = C2 + C3U + Ci{u 2 + 4/i), (j) = 2csh + Ac^uh, 




(5) 
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The maximal Lie invariance algebra g of system (1) with H ^ is also infinite dimensional and 
spanned by the vector fields 

V = td t + xd x , G = td x + d u , C(f,g) = f(u,h)d t + g(u,h)d x . (6) 

The first operator is associated with scaling symmetries in t and x and the second operator 
gives rise to Galilean boosts. Similarly to the usual one-dimensional shallow-water equations, 
the operator C(f,g) indicates that system (1) can be linearized to system (5) by the hodograph 
transformation in which f = t and g = x are the new unknown functions of the new independent 
variables u and h. 

We have verified our computations, leading to the maximal Lie invariance algebra g, using 
the package DESOLV for finding Lie symmetries [35]. 

Besides Lie symmetries, also discrete point symmetries of systems of differential equations 
can be useful, especially in the classification of optimal lists of inequivalent subalgebras or the 
construction of structure-preserving low-dimensional dynamical systems, see e.g. [5,8]. Discrete 
point symmetries can be computed systematically using the direct method or, often more con- 
venient, a version of the algebraic method proposed in [18] and refined in [6, 13] for systems 
admitting infinite dimensional maximal Lie invariance algebras. In the present case, it is found 
that system (1) admits two independent discrete point symmetries (up to composition with 
continuous symmetries and each other), which are given by the pairwise alternation of signs in 
t, x and x, u, respectively. 

Remark 1. In [34] the authors stated that the maximal Lie invariance algebra of system (1) is 
spanned by the four vector fields 

d t , d x , td x + d u , td t + xd x . 

Comparing this algebra, which we denote by g 1 , with the algebra spanned the the vector 
fields (6), it is obvious that g 1 is only a subalgebra of g since the first two vector fields are special 
instances of the parameterized vector field £(f,g), for which (/,<?) = (1,0) and (/, g) = (0, 1), 
respectively. It is thus also missed in [34] that system (1) can be linearized using a hodograph 
transformation. 

In the following, we denote by G 1 the subgroup of the maximal Lie invariance pseudogroup 
G of the modified shallow- water equations (1) that is associated with the subalgebra q . 

3 Classification of subalgebras 

The computation of exact solutions of system (1) was of central importance in [34]. This was 
done with the method of group-invariant (or Lie) reduction. These reductions were constructed 
without a classification of optimal lists of inequivalent subalgebras. Unfortunately such classifi- 
cations are commonly omitted in the literature without justification. This topic is extensively 
discussed in e.g. [26,27], where it is pointed out that only group-invariant solutions stemming 
from inequivalent subalgebras of the maximal Lie invariance algebra of a system of differential 
equations are guaranteed to be inequivalent. This means that they cannot be related to each 
other by means of a symmetry transformation. 

Hence, the first step in the construction of exact solutions of a system of differential equations 
using Lie reduction should be the classification of subalgebras of appropriate dimensions. In the 
present case, system (1) is a system in (1 + 1) dimensions, so it is sufficient to carry out the Lie 
reductions with respect to one-dimensional subalgebras of g, which then yield systems of ordinary 
differential equations. In the following, we classify both subalgebras of the infinite dimensional 
maximal Lie invariance algebra g and the subalgebra g 1 found in [34], as no classification of 
one-dimensional subalgebras of g 1 was given in [34]. 
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The commutation relations between elements of g are 



[V,g}=0, [C(f,g),V]=C(f,g), [Q, £(/, g)] = C(f u , g u - f), [£(/i, 9l ), £(/ 2 , g 2 )] = 0, 

where the pairs (f,g), (/i,<?i) and (72,52) are solutions of system (5). The nonidentical basis 
adjoint actions of the maximal Lie invariance pseudogroup G on the generating vectors fields V, 
Q and C(f,g) then are 

Ad(e e2? )£(/, 9 ) = e s C(f,g), Ad(e^) )v = V _ eC (f, g), 
Ad(e eG )£(/, <?) = £(/',</), Ad{e^^) )g = g + £ £(f u ,g u - f), 

where /' = f(u — e, h) and g' = {g{u — e, h) + ef(u — e, h)). For each pair of generating vector 
fields (v, wo) of g, these adjoint actions can be computed either through the Lie series [26], 

00 n 

w(e) = Ad(e £V )w := V ^-(ad v)"w , 

n=0 

or through a pushforward of the vector fields from g by symmetry transformations from the 
maximal Lie invariance pseudogroup G, see [4, 12] for more details. The second method is 
particularly suitable for computations involving infinite dimensional Lie algebras. 

Now that we have the nonidentical adjoint actions at hand, we can proceed with the classifi- 
cation of the one-dimensional subalgebras of q. To accomplish this classification, we start with 
the most general element of a one-dimensional subalgebra from g, 

v = aV + bg + £(f,g), 

where a, b G M and (/, g) is an arbitrary but fixed solution of the system (5). We subsequently 
simplify v as much as possible by applying the adjoint actions that we found previously. As the 
classification is rather short, we will give it explicitly here. 

(i) If a ^ 0, we can scale a = 1 and use the adjoint action Ad(e e ^'^) to cancel C(f,g) 
from v. No further simplifications are possible and so the simplified form of v is v = T> + bQ. 

(ii) For a = and b ^ 0, we can set 6 = 1 and the adjoint action Ad(e £C ^'^) can be used 
again to cancel C(f,g), and thus v = Q in this case. 

(iii) The final one-dimensional subalgebra is spanned by v = C(f,g). We note that algebras 
spanned by (C(f,g)) and (C(f,g)} are equivalent if the pairs of functions (f,g) and (f,g) differ 
by a constant multiplier or a shift of u. 

Collecting the results obtained, we have thus proved the following theorem. 

Theorem 1. An optimal list of inequivalent one- dimensional subalgebras of the maximal Lie 
invariance algebra q of the modified shallow-water equations (1) consists of the subalgebras 

{V + aG), (G), (J0(f,g)), 

where a£l and (f,g) is an arbitrary but fixed solution of the system (5). 

In principle, Theorem 1 tells us which subalgebras should be used for the Lie reductions of 
system (1). At the same time it is useful to classify the subalgebra g 1 found in [34], for which 
no optimal system of one-dimensional subalgebras was given. This is a simple task, since the 
problem of classifying subalgebras of real four-dimensional Lie algebras is already completely 
solved [28]. See also [25,32] for classifications of real algebras of dimensions up to six. For Lie 
algebras of dimension higher than six exhaustive classification results exist only provided that 
the Lie algebra possesses a special structure: see the review in [32] for more details. 

In the present case, the algebra g 1 is isomorphic to the algebra A® 8 in the classification [32], 
the nilradical of which is (dt,d x ,g), which is isomorphic to the three-dimensional Heisenberg 
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algebra. The one-dimensional subalgebras of the algebra A® 8 can be presented as follows in a 
way suitable for Lie reduction: 



(v + a g), (g + sd t ), (d t + sd x ), (d x ), (7) 

where a £ K and 5 G {0, 1}. It is important to note that we cannot set 5 = in the second 
algebra of the above list as the internal equivalence in g 1 is weaker than the equivalence in g. 



4 Exact solutions of the modified shallow-water equations 

We now present the associated Lie reductions obtained using the optimal list of one-dimensional 
subalgebras of g given in Theorem 1. We solve the arising systems of ordinary differential 
equations whenever possible. 

(%) Subalgebra (T> + aQ). The reduction ansatz in this case is u = u{p) + a hit, h = h(p), 
where p = x/t — a hit + a is the new independent variable. Plugging this ansatz into system (1) 
reduces the two equations to the system of nonlinear ordinary differential equations 



H 



a - pit' + uv! + ( 1 + ^- ) h' = 0, -ph! + uh' + hv! = 0, 



where here and in the following a prime denotes the derivative with respect to p. As this equation 
is still too complicated to be solved in closed form, one could resort to numerical integration for 
finding a solution of the above nonlinear system and then extending the numerical solution to a 
solution of the original modified shallow-water equations. 

(%%) Subalgebra (Q). The ansatz for Lie reduction in this case is u = u{p) + x/t, h = h{p), 
where p = t. The system of modified shallow-water equations (1) then reduces to 

, u ~. h 

u' + _ = 0, h' + - = Q. 
p p 

The solution of this system is u = c\/p and h = o±jt, where c\, C2 E R, giving rise to the Galilean 
invariant solution 

x + ci c 2 

« = — , h = J (8) 

of the modified shallow- water equations (1). 

(Hi) Subalgebra (C(f,g)). For this algebra a suitable reduction ansatz is u = u{p), h = h(p), 
where p = fx — gt. Implicitly differentiating this ansatz with respect to t and x and solving the 
resulting algebraic system for the required derivatives ut, u x , h t and h x yields 

gu P fu p gh p flip 
u t = =r, u x = ——, h t = h x - 



D ' D ' D ' D ' 

where the precise expression for D is not required. Plugging this ansatz into the modified 
shallow-water equations leads to 

-QUp + fuup + / ^1 + hp = 0, -gh p + fuh p + fhu p = 0, 

where / = f(u, h) and g = g(u, h). The reduction with respect to the algebra (C(f,g)) includes 
several physically relevant solution ansatzes. For example, in the case of / = 1 and g = const 
(which obviously is a solution of the system (5) as required), the reduction ansatz coincides 
with that of a traveling wave solution. Similarly, for (/, g) = (1,0) (or (/, g) = (0, 1)) we obtain 
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the stationary (or space independent) reduced system. In the present case, all these reduced 
equations admit only the constant solution u = c\ and h = C2 and thus are not of great interest. 
Less trivial but implicit solutions can be found by choosing / and g in a way such that they 
explicitly depend on u and h. 

We now review succinctly the Lie reductions and results obtained in [34]. As noted in the 
previous section, no optimal list of one-dimensional subalgebras of their algebra g 1 was used 
to systematically carry out reductions. This gave rise to overly complicated reduction ansatzes 
in [34]. Moreover, these ansatzes could lead to solutions which in fact can be mapped to each 
other using a symmetry transformation, i.e. that are equivalent. 

Three Lie reductions were presented in [34]. The first reduction (Case la) is carried out with 
respect to the algebra (a\D + a^Q + a^d x + a±dt), where a± ^ 0. In view of our classification 
results and the optimal system (7), it is obvious that this case is equivalent to our case (i), i.e. 
the algebra (T> + aQ). Thus, the ansatz found in [34] should reduce to the ansatz presented 
above for case (i) when a\ = 1 and 03 = 04 = 0. However, this is not the case. In fact, it can 
be directly checked that the ansatz in [34] for this case is incorrect. The correct (but overly 
complicated) ansatz for the algebra (a\D + CL2G + a^d x + a±dt) would be 



The second reduction (Case lb) rests on the algebra (a2<5 + a^,d x + a^dt). It is equivalent to 
our case (ii), i.e. it is needless to assume 03 7^ 0, 04 7^ 0. However, it is crucial to point out 
that one cannot put 04 = using only the weaker equivalence relation imposed by the adjoint 
action associated with the subgroup G 1 of G on the subalgebra g . This is why it is instructive 
to compare the second case in the optimal list of subalgebras of the maximal Lie invariance 
algebra g given in Theorem 1 with the second case in the optimal list of subalgebras of the 
subalgebra g 1 in (7) to each other. On the other hand, even with the adjoint action of G 1 on 
g 1 , it would be possible to set 03 = 0. 

The last reduction (Case II) employs the algebra {a,2G+a$d x ). The associated group-invariant 
solution is equivalent to our case (ii) i.e. it is again needless to assume 03 7^ 0, and the exact 
solution found in [34] can be obtained from the Galilean invariant solution (8) through re-scaling 
of t and x and shifting of t. 

To sum up, the three reductions carried out in [34] are either incorrect (Case la) or equivalent 
to a solution that could be obtained from a simpler reduction ansatz (Case lb and Case II). From 
the more general point of view, if one would seek only reductions based on the subalgebra g 1 , 
the optimal list (7) should be used to find the proper reduction ansatzes. 

Remark 2. Although the general strategy for the construction of invariant solutions is to use 
inequivalent subalgebras to derive inequivalent solutions, the procedure of generating new solu- 
tions from known ones for system (1) could yield interesting results as well since this procedure 
for system (1) can involve the hodograph transformation. This is why even the use of equivalent 
subalgebras independently will make sense for this system if one is able to construct explicit so- 
lutions and avoiding the use of the hodograph transformations in finding these solutions. On the 
other hand, in [34] this construction would be trivial as the transformations considered in [34] 
are project able. 

So far, we have not used the property that (1) is linearized to (5) by the hodograph trans- 
formation. We now investigate the possibility of finding exact solutions of (5) that can be used 
to obtain exact solutions of the initial system of modified shallow-water equations (1). To this 
end, we start by combining the two equations of the system (5) to a single equation for / by 
excluding g, which reads 





(9) 
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A solution for this equation can be found from the separation ansatz / = fi(u) /2(h) , which 
yields the following form for /, 




/ = ^( cisin (v?") +C2COS (v^ 

where c > 0, c\, . . . ,04 G R, and Y& are the Bessel functions of the first and second kind, 
respectively, and b = y/1 — 4c. As it is quite intricate to obtain the corresponding solution for 
g from the system (5) using the above solution for /, we restrict ourselves to finding particular 
solutions of system (5) following from simpler forms of /. 

To give one example, let us set c = 3/16, 02 = 04 = 0. The above solution then reduces to 




/ = 

and the solution for g from the system (5) becomes 




+ c 5 , 

where C5 6 R. This solution cannot be solved for u and h in terms of / = t and g = x globally, 
but could be solved (e.g. numerically) pointwise and would then yield a solution manifold for 
the original modified shallow-water equations. This solution is obviously inequivalent to the 
group-invariant solutions found above. 

Another, somewhat simpler solution of Eq. (9) is 

u 1 

J = ClT + C 2 T + C 3 1i + C 4 

h h 

which, when substituted back into (5) yields the following solution for g, 
g = a ^j- + ^ - \nh\ + c 2 ^ + c 3 Qu 2 - Hlnh - h^j + c 5 

where again c±, . . . ,05 S R. In the special case when 01 = 03 = 04 = 05 = 0, this solution becomes 
the Galilean invariant solution (8) of the modified shallow-water equations found above. If all 
the constants c\, . . . , C5 are non-zero, it is again impossible to solve the above solutions for u 
and h globally in terms of / = t and g = x. 

In the same way one could proceed to find particular solutions of the system (5) and then 
extend them to solution of the original modified shallow-water equations. We will not pursue 
this idea here further though. 



5 Conclusion 

In the present paper we have reconsidered the problem of studying the Lie symmetries of a system 
of modified shallow- water equations (1) recently investigated in [34]. We have derived that the 
maximal Lie invariance algebra of this system is in fact infinite dimensional, containing the 
algebra found in [34] as a subalgebra. The infinite dimensional part of the maximal Lie invariance 
algebra (6) indicates that the system (1) can be linearized using a hodograph transformation, 
the result of which is the system (5). This is not unexpected, as any homogenous first-order 
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system of partial differential equations with two dependent and two independent variables can 
be linearized by a hodograph transformation, provided that the coefficients in the system only 
depend on the unknown functions. As in [34] only a finite dimensional symmetry algebra is 
found, the indication, provided by the infinite dimensional maximal Lie invariance algebra, of 
the existence of a linearizing transformation is missed. 

When it comes to determining exact solutions of the system (1) we observe some common 
errors in finding exact solutions of partial differential equations in [34] . These errors are discussed 
in [21,33]. The most crucial point missed is that the algebra g 1 is not maximal as a symmetry 
algebra of the system (1). Moreover, no optimal list of inequivalent one-dimensional subalgebras 
of the algebra q 1 was constructed. As the algebra q 1 is only four-dimensional this would have 
been a simple task as for four-dimensional Lie algebras the problem of classifying all inequivalent 
subalgebras is already completely solved. If one omits the construction of such an optimal list 
of inequivalent subalgebras prior to carrying out the Lie reductions, then one might end up with 
unnecessarily complicated reduction ansatzes. This then leads to needlessly complicated exact 
solutions for the reduced systems. We observed this with the exact solutions found in [34]. In 
another light, carrying out Lie reductions without reference to the optimal list of inequivalent 
subalgebras can lead to equivalent exact solutions, i.e. solutions that are related through a 
symmetry transformation. 

Finally, as the linearizing hodograph transformation was missed in [34] the authors could not 
construct exact solutions of the linearized modified shallow- water equations (5), some of which 
can be mapped to exact (and non-Lie) solutions of the modified shallow-water equations (1). 
This is another possible source for finding exact solutions of system (1). 
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